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ABSTRACT 

We present position-velocity strip maps of the Galactic Center region in the 
CO J = 7 ^ 6 and J — 4 ^ 3 transitions observed with the Antarctic Sub- 
millimeter Telescope and Remote Observatory (AST/RO) located at Amundsen- 
Scott South Pole Station. Emission from the two rotational transitions of ^^CO 
was mapped at 6 = 0° for 3.5° > i > —1.5°, on a 1' grid with a FWHM beamsize 
of 58" at 806 GHz and 105" at 461 GHz. Previous observations of CO J = 4 ^ 3 
(Martin et al., in preparation) and [CI] (Ojha et al. 2001) emission from this re- 
gion show that these lines are distributed in a manner similar to CO J = 1 — >■ 
(Stark et al. 1987); the (CO J = 4 ^ 3)/(C0 J - 1 ^ 0) line ratio map 
is almost featureless across the entire Galactic Center region. In contrast, the 
CO J = 7 — >^ 6 emission from the Galactic Center is strongly peaked toward 
the Sgr A and Sgr B molecular complexes. A Large Velocity Gradient (LVG) 
analysis shows that aside from the two special regions Sgr A and Sgr B, the 
photon-dominated regions within a few hundred parsecs of the Galactic Center 
are remarkably uniform in mean density and kinetic temperature at n = 2500 
to 4000 cm-3 and T = 30 to 45 K. The (CO J = 7 ^ 6)/(C0 J = 4 ^ 3) line 
temperature ratios near Sgr B are a factor of two higher than those observed 
in the nuclear region of the starburst galaxy M82 (Mao et al. 2000), while the 
CO (J = 7 — 6)/C0(J = 4 — 3) line temperature ratios around Sgr A are 
similar to M82. The line ratio on large scales from the Galactic Center region is 
an order of magnitude less than that from M82. 

Subject headings: Galaxy: abundances — Galaxy: center — ISM: general — 
ISM: molecules 
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1. Introduction 

Observations of photons emitted by the various rotational transitions of the ground 
vibrational state of carbon monoxide (CO) are the primary means of studying molecular gas 
in the Galaxy. These spectral lines occur at frequencies of J x 115 GHz, for the transition 
from the J to the J — 1 rotational state. Numerous galactic surveys (Combes 1991) have 
studied the lowest-frequency J = 1 — >■ hne. The brightness of this spectral hne is roughly 
proportional to the total molecular column density within the telescope beam (Liszt 1984), 
provided that the molecular gas has a relatively low density (n < 3 x 10^ cm~^) and column 
density {N < lO^^cm"^). Surveys of the J — 1 ^ line therefore give an indication of the 
extent and distribution of molecular material. 

A more complete picture of the thermodynamic state of the molecular gas can be 
gained through observations of other, higher- J transitions of CO and the various rotational 
transitions of the isotopically-substituted species ^^CO and C^®0. If the brightness of 
several of these lines is known, models of the excitation and radiative transfer can be used 
to solve for the density, temperature, and coohng rate of the molecular gas. To avoid 
degeneracy in this solution, it is valuable to have observations of a line transition from a 
J-state that is sufficiently high in energy such that it is only weakly populated. CO in 
interstellar molecular gas is typically in a thermodynamic state where the low- J transitions 
are in approximate thermal equilibrium, so that all of the rotational transitions below some 
J-level have roughly the same excitation temperature as the J = 1 ^ transition. The 
observable quantity that is usually compared with theory is the ratio of line brightnesses. If 
only low- J lines are observed, it is found that the ratio of these line brightnesses are only 
weakly dependent on excitation temperature and that maps of all the low- J transitions look 
similar, regardless of variations in excitation temperature across the map. If, however, we 
observe a transition from a high- J level that is not strongly excited, it will be found that 
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the line brightness ratios involving that transition do vary significantly across the map, and 
the value of the excitation temperature can be more readily determined. 

Previous observations with the AST/RO telescope (see Figure 1) have shown that 
the distribution of the CO J = 4 ^ 3 line emission is remarkably similar to that of the 
CO J = 1 — > line throughout the Galactic Center region. Previous observations of CO 
J = 7 — > 6 (Harris et al. 1985) showed that hue to be unusually bright within ±200" of the 
Galactic Center. In this paper, we report observations of the CO J = 7 — > 6 and J = 4 ^ 3 
transitions in a strip map several degrees in extent, and find that the J = 7 — 6 line has 
significantly lower brightness temperature and significantly different spatial distribution 
than the lower- J CO lines. We can therefore produce a new model of the thermodynamic 
state of the CO gas along this 1-dimensional strip. 



2. Observations 

These observations were made from 21 to 25 July 2000 with the Antarctic Submillimeter 
Telescope and Remote Observatory (AST/RO), located at Amundsen-Scott South Pole 
Station. AST/RO is a 1.7 m diameter, offset Gregorian telescope capable of observing at 
wavelengths between 200 iim and 1.3 mm (Stark et al. 2001). 

The facility dual-channel SIS waveguide receiver (Walker et al. 1992; Honingh et al. 
1997) was used in position-switching mode. Emission from the 806.6517 GHz J = 7 — > 6 
transition of ^^CO was mapped at 6 = 0°, 3.5° > / > —1.5° with 1' spacing and a beam 
size of 58". The 461.041 GHz J = 4 — > 3 transition of ^^CO was observed simultaneously, 
with a beam size of 105". The atmosphere-corrected system temperature during the 
observations ranged from 12,000 to 26,000 K at 806 GHz and from 3000 to 4000 K at 
461 GHz. Observing time at each of the 304 observed positions was typically 7 minutes. 
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The intermediate frequency output of each receiver was connected to one of two 
acousto-optical spectrometers (Schieder, Tolls, and Winnewisser 1989) of which 900 x 
1-MHz- wide channels were in use at 806 GHz and 1000 x 1-MHz- wide channels were in 
use at 461 GHz. This configuration provided a velocity range of 225 km at 806 GHz 
and 438 km s"^ at 461 GHz. The channel spacing of 0.67 MHz corresponds to a velocity 
resolution of 0.37 km s~^ at 806 GHz and 0.65 km at 461 GHz. The high frequency 
observations were made with the CO J — 7 ^ 6 line in the lower sideband (LSB). Since 
the intermediate frequency of the AST/RO system is 1.5 GHz, the -^^ Pi hne of 
[C I] at 809.34197 GHz appears in the upper sideband and is superposed on the observed 
LSB spectrum. The local oscillator frequency was chosen so that the nominal line centers 
appear separated by ~100 km s~^ in the double-sideband spectra. Velocity increases with 
spectrometer channel number in opposite directions for the two sidebands; in Figure 4, for 
example, velocity increases to the left for the [CI] line and to the right for the CO. There is 
809 GHz [C I] emission superposed on the 806 GHz CO J — 7 ^ 6 data at most positions 
where the J = 7 — > 6 line is strong. The central channel of the spectrometer was shifted 
for observations at negative longitudes so that the limited bandpass covered most of the 
velocity range of known low- J emission. 

During the time these observations were made, the AST/RO automated cahbration 
system (Stark et al. 1997, 2001) was broken, so that we were not able to make automated 
sky or receiver measurements. We did, however, make manual measurements once a day. 
These show that the receiver gain and noise were stable within a few percent throughout 
the period of the observations. To determine the atmospheric opacity and sky temperature, 
we used data from the NRAO-CMU 350 /im tipper in operation at the South Pole. 
Several times each hour, this instrument makes a calibrated skydip using a broadband 
room-temperature bolometer. During the subsequent observing season, when the automated 
calibration system was working, it was determined that the opacity in the narrow 806 GHz 
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band is related to the opacity measured by the NRAO-CMU tipper through the relation 
(R. A. Chamberlin, private communication): 

TsoeGHz = 1-82 tnrao-cmu — 1-51 . 

The relation used to cahbrate the 461 GHz data was (Chamberhn 2001): 

T461GHZ = 0.55 Tnrao-cmu — 0.25 . 

We used these relations, together with the NRAO-CMU tipper data and thermometer 
measurements of the outside ambient temperature, to calculate the Tgicy value we would 
have gotten if the automated calibration system had been working. These values were 
edited into a copy of the data file, allowing the data to be treated normally by our data 
reduction program, COMB, as described in Stark et al. (2001). 

3. Results and Analysis 
3.1. Strong CO J = 7 — >^ 6 Emission in the Sgr B and Sgr A Complexes 

Longitude-velocity maps of the ^^CO J = 4 — >^ 3 and J = 7 — >■ 6 distributions are 
shown in Figure 2. As can be seen by reference to the (£, b) maps in Figure 1, the strip 
maps at 6 = shown in Figure 2 skirt the outer regions of the strongest CO features, rather 
than passing through their peak positions. Nonetheless, the CO J = 7 — >^ 6 emission is 
concentrated in the Sgr B and Sgr A complexes at velocities between and -1-100 km s~^. 
Emission features at negative velocities in Figure 2b arc contributed by superposed 809 GHz 

— ^■^ Pi [CI] emission from the upper sideband (see §2). 

For comparison, longitude-velocity maps of the ^^CO J = 1 ^ 0, ^"^CO J = 1 — >• 0, 
and CS J = 2 ^ 1 data described by Stark et al. (1987) and Bally et al. (1987, 1988) are 
shown in Figure 3. Again, it is seen that the low-J rotational transitions are broadly and 
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similarly distributed, while the CO J = 7 — > 6 emission is much more spatially confined. A 
direct superposition of the ^^CO J = 4 — > 3 and J = 7 — > 6 maps is shown in Figure 4, 
along with sample spectra at several positions. 

Figure 5 shows maps of ^^CO J = 7 ^ 6 to ^^CO J = 4 ^ 3 {TjljTllr^) and 
^^CO J = 1 ^ to ^^CO J = 1 ^ (T/^o/^iio) line brightness temperature ratios 
as a function of longitude and velocity in the Galactic Center region. The data were all 
convolved to the resolution of the ^^CO J = 1 — > map. The data in this figure fall into 
three categories: (1) regions where the signal-to-noise ratio is high for both the numerator 
and the denominator, (2) regions where the signal-to-noise ratio is low for the numerator 
but high for the denominator, and (3) regions where the signal-to-noise ratio is low for 
both the numerator and denominator. The fourth possibility, where the signal-to-noise 
ratio is high for the numerator but low for the denominator, is avoided by placing the 
higher signal-to-noise data in the denominator. Category 1 data appear as smooth blue, 
green, yellow, orange, or white areas of the map. In these regions, the temperature ratio is 
well-determined, and has a fractional error only a little larger than that of the temperature 
data constituting the numerator. Category 2 data appear as mottled black and blue, since 
the value of the numerator is small and dominated by noise, while the denominator is some 
well-determined number ([zero ± noise] /[number]); the values of the ratio in these regions 
are near zero, with large fractional error. Category 3 data have a "salt and pepper" black 
and white appearance as the ratio varies wildly over all possible values, both positive and 
negative, since both the numerator and denominator have noisy values near zero. Both the 
value of the ratio and its error are large. In the upper panel of Figure 5, category 1 data 
can be found in the vicinity of Sgr A and Sgr B; category 2 data is found between Sgr A 
and Sgr B and to the left of Sgr B; category 3 data occupies the periphery of the figure. 

Typical values of Ty^g/T^^g « 0.3 ± 0.05 at points near Sgr A. At positions on the 
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strip near Sgr B, typical values are higher: T^'^^/Tl^^ Ri 0.6 ± 0.05, over the velocity range 
50 to 85 km . Average hne temperature ratios near Sgr A and Sgr B are summarized 
in Table 1. The value of T^'^^/Tl'^^, i^^ar Sgr B is a factor of two higher than the ratio 
observed in the nuclear region of the starburst galaxy M82 (Mao et al. 2000). Figure 5 
shows that aside from the regions near Sgr A and Sgr B, T^^g/T^^a ~ 0.12 ± 0.08 in the 
region bounded roughly by 0.1 < ^ < 1.5 and —10 km s~^ < v < 110 km s~^. This is the 
area that shows significant CS J = 2 ^ 1 emission in Figure 3 and is the area that Binney 
et al. (1991) identify with X2 orbits in the galactic bulge. The 300 parsec ring is the rough 
parallelogram-shape between (.=—\° and £=1.8°, and velocities from —200 to +200 km s~^. 
The Clump 2 molecular cloud complex is located near £=3°, at velocities from +20 to +150 
km s~^ . Both these regions show T^^^/T^^^ < 0.1. 



3.2. LVG Analysis Using the Line Ratios 

Using the line ratios, we can estimate the kinetic temperature, Tkin, and the number 
density of molecular hydrogen, n(H2), through a large velocity gradient (LVG) radiative 
transfer analysis (Goldreich & Kwan 1974). The LVG approximation simplifies radiative 
transfer by the assumption that an emitted spectral line photon can only be absorbed 
"locally" , within a small region whose velocity is similar to the point of emission. This 
approximation is robust, in the sense that the results of LVG models are often reasonably 
accurate, even when reality violates the assumptions underlying the models (Ossenkopf 
1997). Our LVG radiative transfer code simulates a plane-parallel cloud geometry. It 
uses the CO coUisional rates from Turner (1995) and newly-derived values for the H2 
ortho-to-para ratio (~ 2) and the coUisional quenching rate of CO by H2 impact (Yan, 
Balakrishnan, & Dalgarno, in preparation). The model has two input parameters: the ratio 
of ^^CO to ^^CO abundance, and the ratio X{CO)/W , where X{CO) is the fractional 
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CO abundance parameter and W denotes the velocity gradient. The abundance ratio 
^^CO/^^co is taken to be 25 in the Galactic Center region (Langer & Penzias 1990, 1993). 
We take X{CO)/VV = lO'^-^ pckni'^ s, assuming that the ^^CO/Ha ratio is 10"^ and the 
velocity gradient of Galactic Center clouds is typically 3 to 6 km s~^pc~^ (Dahmen et al. 
1998). 

The hne ratios for the LVG analysis and the physical parameters derived are 
summarized in Table 1. The solution of the radiative transfer equation gives T^^g/T^^g 
and T/^q/T/^^q as a function of the total molecular hydrogen volume density for different 
values of the gas kinetic temperature, as shown in Figure 6. For each observed point, we 
can invert these functions to determine the kinetic temperature and molecular hydrogen 
volume density corresponding to the observed line ratios. 



3.3. Temperature and Density Variations in the Galactic Center 

The LVG model shows the variation of kinetic temperature across the Galactic Center 
region. Sgr B, with Ty^g/T^^g ^0.6, has a kinetic temperature of 72 (±2) K, the highest 
kinetic temperature among the Galactic Center molecular clouds on our strip map. This 
derived kinetic temperature is consistent with the rotational temperature derived from the 
NH3(J,K)=(3,3) ^ (1,1) data (Morris 1989) and ^^CO J =1^0 survey (Oka et al. 1998). 
The kinetic temperature at the position nearest Sgr A is 47 (±3) K, lower than Sgr B. If we 
exclude these two regions, we find that the mean kinetic temperatures is ~35 ±10 K within 
the inner 500 parsccs of the Galaxy. In comparison, the dust temperatures in the inner 200 
pc are ~21±2 K (Pierce-Price et al. 2000). 

The observed line intensities of CO J=l— *0, J = 4^3, and J = 7 ^ 6 lines, 
together with the LVG values for the kinetic temperature, allow us to estimate the column 
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densities of all the J-states of CO. The CO J = 4 — > 3 line is the most luminous CO 
transition from the Sgr A complex, whereas the J = 5 — > 4 line is the most luminous CO 
transition from the Sgr B complex. 

Adopting diameters of 45 pc for Sgr B and 30 pc for Sgr A from inspection of ^^CO 
and CS maps (Bally et al. 1987), and assuming uniform thermodynamic properties within 
those volumes, we can estimate the properties of those clouds as a whole. The total CO 
J = 7 — >■ 6 luminosity integrated over the Sgr B complex is ~ 5 x 10^^ ergs s~^ and 
integrated over the Sgr A complex is ~ 1.5 x 10'^^ ergs s~^. We estimate the total CO 
rotational hne cooling to be Leo ~ 2.9 x 10^^ ergs s~^ from Sgr B and Lqo ~ 2.0 x 10^^ 
ergs s~^ from Sgr A. The CO J = 7 — >■ 6 transition accounts for 7.5% of the total CO 
luminosity from the Sgr A complex (Table 2). At the Sgr B complex, with its higher gas 
temperature, the CO J = 7 — 6 transition contributes more than 17% of the total CO 
luminosity. The total estimated energy released by the CO rotational lines toward the 
Sgr B cloud is, however, similar to that of Sgr A. 

Using the derived molecular hydrogen volume density from our LVG analysis and 
assuming spherical geometry, we can estimate the masses of the Sgr B and Sgr A molecular 
clouds. Multiplying the LVG-calculatcd density by the size of Sgr B, we find its mass is 
approximately 10^ M©, consistent with estimates of ~8 x 10^ Mq, derived using dust 
continuum emission (Lis, Carlstrom, & Keene 1991; Gordon et al. 1993). The CO 
J = 7 — >■ 6 emitting region towards the Sgr A complex has a mass of ~ 0.4 x lO^M©. 

4. Summary &; Conclusions 

Wc have presented maps of the Galactic Center region in the CO J = 7 — > 6 and 
J = 4 — 3 transitions observed with the Antarctic Submillimeter Telescope and Remote 
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Observatory (AST/RO) located at Amundsen- Scott South Pole Station. Comparing the 
CO J = 4 ^ 3, CO J = 1 ^ 0, and ^^CO J = 1 ^ maps reveals that the CO J = 7 ^ 6 
emission is concentrated toward the Sgr A and Sgr B complexes. Using an LVG model, we 
find that the gas kinetic temperatures of the Sgr B and Sgr A complexes are 72 K and 47 
K and the gas densities of those regions are 10^"^^ cm~^ and 10^'^ cm~^. Excluding those 
regions, we estimate the mean kinetic temperature and density throughout the Galactic 
Center are 30 to 45 K and lo^ ^io.os ^^-3 rpj^g qq j — ^ ^ q transition accounts for 7.5% 
of the total CO luminosity from the Sgr A complex. At the Sgr B complex, with its higher 
gas temperature, the CO J = 7 — > 6 transition contributes more than 17% of the total CO 
luminosity. 

The tragic death of Rodney Marks, who was the year 2000 AST/RO Winterover 
Scientist, occurred just before this observational program. The remaining winterover crew 
from the Center for Astrophysical Research in Antarctica — Gene Davidson, Greg Griffin, 
David Pernic, and John Yamasaki — continued AST/RO operations in tribute to Rodney's 
memory, allowing these observations to be made. The AST/RO group is grateful for the 
logistical support of the National Science Foundation (NSF), Antarctic Support Associates, 
Raytheon Polar Services Company, and the Center for Astrophysical Research in Antarctica 
during our polar expeditions. SK thanks Min Yan for LVG analysis and Wilfred Walsh for 
his helpful comments on the manuscript. This work was supported in part by United States 
National Science Foundation grant DPP88- 18384, and by the Center for Astrophysical 
Research in Antarctica and the NSF under Cooperative Agreement OPP89-20223. 
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Table 1: Parameters and Results of LVG Model Calculations 

Parameter Sgr A Sgr B 

T^tjTlt^i 0.3 ± 0.05 0.6 ± 0.05 

Tl^jTl^Q 5.7 ± 1.6 5.0 ± 1.0 

Velocity Range (km s"^ ) 10 - 80 50 - 85 

n(H2) (cm-3) iQZ.^±om io4.i5±o.o5 

Tk^n (K) 47 ± 3 72 ± 2 

Mass H2 {Mq) 0.4 X 10^ 1.0 xlO^ 

Size (pc) 30 45 



Table 2: CO & CI Luminosity 

Sgr A Sgr B 

L(CO 7^6) (ergs s-i) 1.50x10^6 5.07x10^6 

L(CO 4^3) (ergs s-i) 4.73x10^6 5.59x10^6 

L(CO total) (ergs s-i) 2.00x10=^^ 2.92x10^^ 

L(C I total) (ergs s-i) 5.46x10=^^ 8.96x10^5 
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J? (degrees) 

Fig. 1. — Top: Integrated intensity map of the ^^CO J = 4 3 emission observed with the 
AST/RO telescope (Martin et al., in preparation). Bottom: Integrated intensity map of the 
^^CO J = 1 ^ emission from the same region (Bell Labs 7m data, Stark et al. 1987). Ten 
equally spaced contours are drawn at the levels shown in the colorbar above. 
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Fig. 2. — Longitude- velocity maps of CO emission observed toward the Galactic Center 
with AST/RO. Top: J = 4 ^ 3. Bottom: J = 7 ^ 6. The emission components nominally 
appearing near —50 km s^^ in the bottom map are from the 809 GHz [CI] line (see text). 
The Sgr A cloud is near i = 0° and the Sgr B cloud is near i = 0.67°. The dark horizontal 
bands in the top map near f = to —60 km s~^ are due to absorption in foreground gas. 
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Galactic Longitude (degnees) 

Fig. 3. — Longitude- velocity maps of ^^CO (J = 1 — > 0) emission [top), from Bell Labs 7m 
data described by Stark et al. (1987); ^^CO (J = 1 ^ 0) (middle) and CS (J = 2 ^ 1) 
emission [bottom), from data presented by Bally et al. (1987). All three data sets are shown 
smoothed to 3' resolution. The ^^CO data for £ > 0.9° and i < -0.5° were observed with a 
spacing of 6', leading to a coarse appearance in these regions. 
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Fig. 4. — Sample ^^CO J = 7 ^ 6 {blue lines) and ^^CO J = 4 — > 3 {red lines) spectra 
observed toward 4 different positions, as marked. Center map shows the ^^CO J = 7 ^ 6 
contours (in blue) superposed on the ^^CO J = 4 3 longitude- velocity map (in red). 
Superposed emission from the 809 GHz [CI] line in the upper sideband is shown in light 
blue. The ^^CO J = 7 ^ 6 spectra are shown after Hanning smoothing. The spectra near 
Sgr B {i = 0.67°) have higher signal-to-noise ratio due to significantly greater integration 
time than the other, more typical spectra. 
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Fig. 5. — Longitude- velocity maps of the ratios of line temperatures of CO transitions in the 
Galactic Center region: (a) Tj^q/T^^^ (this paper), (b) T^^q/T^'^q (data from Stark et al. 
1987). 
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Fig. 6.— Large Velocity Gradient (LVG) model results for X{CO)/W = lO-^-^pckm-^ s 
{X{CO) is the fractional CO abundance parameter). Solid and dashed lines show constant 
values of T^^^/Tl^^ and of T^^q/T^^q, as marked. EUiptical symbols represent the data with 
one sigma errors, from Table 1. The analysis uses newly derived values for the H2 ortho-to- 
para ratio and for the collisional quenching rate of CO by H2 impact (Yan, Balakrishnan, & 
Dalgarno, in preparation). 



